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Abstract. We discuss the rheology experiments on nematic elastomers by Martinoty et al. in the light of 
theoretical models for the long-wavelength low-frequency dynamics of these materials. We review these 
theories and discuss how they can be modified to provide a phenomenological description of the non- 
hydrodynamic frequency regime probed in the experiments. Moreover, we review the concepts of soft and 
semi-soft elasticity and comment on their implications for the experiments. 
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Nematic elastomers are crosslinked rubbers with the 
uniaxial symmetry of a nematic liquid crystal fj]. One 
could argue, therefore, that they are simply uniaxial rub- 
bers. A nematic elastomer can, however, be distinguished 
from a simple uniaxial elastic medium by the fact that, at 
least in an idealized limit, it can form via a spontaneous 
phase transition from an isotropic phase. Since this phase 
transition breaks the continuous rotational symmetry of 
the isotropic phase, it has an associated Goldstone mode 
whose manifestation is the vanishing of the elastic mod- 
ulus C5 measuring the elastic energy of strains in planes 
containing the anisotropy axis |2HH]- Thus, the ideal ne- 
matic elastomer has a "soft" elasticity compared to a tra- 
ditional uniaxial elastic medium in which C5 is nonzero. 
In practice, ideal monodomain nematic elastomers do not 
form in the absence of some aligning field, which is usually 
produced by weakly crosslinking a sample in the isotropic 
phase, stretching to produce a uniaxial configuration, and 
then crosslinking again [I] . The resultant material is weakly 
anisotropic in the high-temperature paranematic phase 
and more anisotropic in the low-temperature nematic phase. 
Since the low-temperature phase arises from an already 
aligned phase and does not break rotational symmetry, it 
does not exhibit ideal soft elasticity. Rather, it exhibits 
"semi-soft" elasticity in which C5 is small but nonzero. 
Thus, for small strains, a semi-soft nematic elastomer is 
truly a uniaxial solid. At higher strains, however, it ex- 
hibits properties, such as a nearly constant stress for in- 
creasing strain, characteristic of a soft nematic elastomer 

u 

The static elastic properties and phase transitions of 
both soft and semi-soft nematic elastomers are reason- 
ably well understood |H2l6l7ftj| . provided the effects of 
random stresses 2,9 can be ignored. They can be de- 
scribed either in terms of models involving strain only or 
in terms of models with coupling between a traditional 



Maicr-Saupe-de-Gcnnes nematic symmetric-traceless or- 
der parameter Qy and strain. Considerable progress has 
also been made toward understanding the dynamics of 
nematic elastomers [10111 1U12| again in the limit in which 
random stresses and fields can be ignored. As in a stan- 
dard elastic medium, purely hydrodynamical equations, 
which describe all modes with frequencies u> smaller than 
the smallest characteristic inverse decay time r _1 when 
wavenumber q tends to zero, involve only the displace- 
ment field u and not the nematic director n, which re- 
laxes to the local stain in a nonhydrodynamic time r n . 
Even though the nematic director is not strictly speak- 
ing a hydrodynamic variable in either soft or semi-soft 
nematic elastomers, it is of some interest to develop phe- 
nomenological equations describing both nematic director 
and elastic displacement. Two distinct derivations [TTlllI 
IT2"] lead to identical sets of coupled director-displacement 
equations, which we will refer to as the NED (nematic- 
elastomer dynamic) equations, that reduce to the rigorous 
hydrodynamic equations in terms of displacement only at 
frequencies Lor n <C 1 and to the standard equations of 
nematohydrodynamic when elastic rigidity is turned off. 
Rheological measurements of the complex modulus 



G*{uj) = G'{lu)-iG"{uj) 



(1) 



provide a useful experimental probe of the dynamical prop- 
erties of viscoelastic systems. They typically probe a wide 
range of frequencies that stretch from the low-frequency 
hydrodynamic regime to well above it. Thus, interpreting 
these experiments requires theoretical models of dynamics 
beyond hydrodynamics. A nematic elastomer is character- 
ized in general by relaxation times associated with director 
relaxation and with other modes, which we will simply re- 
fer to as elastomer modes. In their simplest form, the NED 
equations assume the director relaxation time r n is the 
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longest non-hydrodynamic relaxation time, and that it is 
well separated from the longest elastomer time te- They 
describe dynamics for both toT n <C 1 and u)T n > 1 and 
provide non-trivial predictions about G*(w), in particular 
the appearance of a plateau in G'(w) and an associated 
dip in G"{uj) for r" 1 < to < t e y . Experiments [T11TH11TCI 
I17U18| to date do not exhibit the plateau predicted by the 
NED equations. This has led Martinoty et al. to question 
their validity. 

Here we will review both the hydrodynamical and NED 
equations for nematic elastomers and discuss how they can 
be modified to provide a phenomenological description of 
the non-hydrodynamic frequency regime probed in rheol- 
ogy experiments. Moreover, we review the concepts of soft 
and semi-soft elasticity and comment on their implications 
for rheology experiments. 

We begin with a review of the dynamics of traditional 
isotropic solids, of which rubber is a particular example. 
The elastic free energy of an isotropic solid is character- 
ized by a bulk modulus B and a shear modulus /i, and 
dissipation is characterized by a bulk viscosity £ and a 
shear viscosity r/, which are both frequency independent 
at low frequencies. The complex shear modulus in the hy- 
drodynamic limit is simply 

G*(lo) = (JL-iur). (2) 

This form is rigorous in the hydrodynamic limit where 
lot <C 1. /i and rj can vary widely in magnitude from 
system to system. Nothing more can be said about the 
behavior of G* (w) without some more microscopic model 
of the dynamics of a particular system. 

The Rouse model provides a simple description of 
the dynamics of polymeric melts. The dynamics of even 
isotropic elastomers and rubbers is far more complex than 
that of melts, and the Rouse model does not necessar- 
ily provide a good physical description of their dynami- 
cal modes. Nevertheless, experimental curves for G*(w) in 
nematic elastomers are very similar to those predicted by 
the Rouse model, showing in particular a region in which 
G'{lo) and G"(lo) track each other with a w° behavior with 
a « 1/2. It is thus legitimate to use Rouse- model pre- 
dictions as fitting curves for nematic elastomers, whether 
or not it provides a microscopically correct description 
of their dynamics. This is the point of view we will take 
here. We parameterize the complex storage modulus in 
the isotropic phase of nematic elastomers as 

G*(lo) = G [l + ah R (-iujT E )], (3) 

where Go = n x kT, with n x being the number density of 
polymer strands between crosslinks, and a — Yl p =iP~ 2 - 
hn{x) = x Jr(x) is a Rouse function with 

fn{x) = — jj • (4) 

In the limit N — > oo, the summations in Eq. can be 
evaluated via the Poisson summation formula with the 
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h R (x) = , =- {7rV^cosh(7r v / a;) - smYi(ir^fx) \ . 

7r z siaxiyK^/x) 

(5) 

When lute <C 1, this expression reduces to Eq. J3J) with 
/i = n x kT, and r\ = (tt 2 /6)hte- When t^ 1 < uj < r^ 1 , 
G'{lo) ~ G"(w) ~ V^- 

A uniaxial solid is characterized by an elastic energy 

f ( G G 
Hu= I d ^ x \~^ + c * u zz uu + — u% 

+ C 4 ul b + C«ul z y (6) 

with five elastic constants, where Uy with i,j =x,y,z is 
the strain tensor and a,b = x, y. The superscript R on 
the elastic modulus for shears in plane containing the 
z-axis indicates that it has been renormalized by director 
fluctuations. In soft nematic elastomers, G|* is zero, and it 
is necessary to add bending terms proportional to (d 2 u a ) 2 
to ensure stability in all directions. In practice all single 
domain nematic elastomers are semi-soft with G,f > 0, 
and the bending terms can be neglected. Associated with 
each of the five elastic constants is a viscosity -qk- The 
complex storage moduli for shears in the planes containing 
the z-axis and in the plane perpendicular to the z-axis are, 
respectively, 

qH = G«-z^, (?) 

G* ± (u) = Ci - iujrji. (8) 

As in isotropic rubbers, we can say nothing more about the 
behavior of Gm(cj) and G* ± (u>) without some microscopic 
model of the dynamics. Since there are surely elastomer 
modes in the anisotropic phase of nematic elastomers, one 
could argue that they are dominant and that Gn and G*j_ 
have exactly the same form as G in an isotropic system, 
Eq. (0), with possibly different values of Go, a and te- 
Indeed, the data presented in Martinoty, et al. appear to 
fit this model with t e ~ 10 -2 , considerably longer than 
the Rouse time in traditional isotropic rubbers. Certainly, 
however, there are extra director modes in nematic elas- 
tomers that are not present in isotropic ones, and a com- 
plete model of dynamics beyond hydrodynamics must in- 
clude them. 

A phenomenological approach to the dynamics of ne- 
matic elastomers including the nematic director begins 
with an elastic free energy that includes both strain and 
the director. Following de Gennes [201 ^ ^ ne energy Ti can 
be written as a sum of the Frank free energy 7i n for the 
director, the strain and strain-director coupling energy 

/( G G 
d 3 x <—u 2 zz + C 2 u zz uu + — u% 

+ C 4 u 2 ab + C 5 u 2 az + ^-Q 2 a + D 2 u za Q a \, (9) 
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where 



Q a = 5n a - - (d z u a - d a u z ) . 



(10) 



When Q a relaxes to its equilibrium value — {D 2 /Dx)u az 
in the presence of strain u az , the energy 7i u . n reduces to 
H„ with C* = C* 5 - D 2 2 /{2D 1 ). 

The phenomenological equations for n and u, which 
can be derived using standard Poisson-bracket approaches 
EU, are [12] 

hi = Xi 3 k dju k - r ^ , (11a) 

pili = X kji dj^- - |^ + uukt djdiiik , (lib) 
dn k Out 

where Vijki is an unrenormalized uniaxial viscosity ten- 
sor with five independent components, r is a dissipative 
coefficient with dimensions of an inverse viscosity, and 

hjk = ^ (Sjj n k + 5j k n k ) + ^ (Sfj n k - 5f k n k ) , (12) 

where 6fj = 6ij — riinj. The complex modulus GfJco) — 
C^(uj) predicted by these equations is 

W»)=C,-i»L + £)-£±p^, (13) 



2T 



2 D\ 1 — iuTi 



where t x = l/(-TDi) and t 2 = -X/(rD 2 ). This result is 
equivalent to the one obtained by Terentjev, Warner and 
coworkers [Tllll| . Equation (|13Jl is a valid expression for 
the complex modulus provided n = r„ is much longer 
than any other nonhydrodynamic relaxation time, includ- 
ing elastomer times, of the system. As can be seen in 
Fig. I^a), it predicts a plateau in G'(uj) for w > t~ 1 . 

In standard polymeric systems, the Rouse time, tr, is 
of order 10 _5 s, and light scattering measurements on ne- 
matic elastomers yield r n ~ 10 -2 [221 - It could, therefore, 
be argued that te ~ t r and that that r n ^> te- However, 
liquid crystalline polymers are more viscous than stan- 
dard polymers, and it is not unreasonable for te to be 
considerably less than 10 _5 s. If t„ and te are of the same 
magnitude, or if r n < te, then the viscosities and decay 
times T\ and t 2 in Eq. I|13fl must be treated as having a 
frequency dependence when ujte is of order 1 or greater. 

Since high-frequency dynamics of the nematic phase 
are expected to be similar to that of the nearly isotropic 
paranematic phase, should have the same behavior as 
the viscosity in the latter phase: 



v b (w) = V5fn{-iuT E ). 



(14) 



In addition, the elastomer modes will affect the relaxation 
of director modes. Indeed, the light scattering experiments 
of Ref. [221 indicate a very broad distribution of director 
relaxation times. There is no reason why many closely 
spaces modes should not lead to a behavior of _T _1 with 
a Rouse-like form, 



ilute) 



(15) 



Cj io 1 



10" 



(a) r„/r £ = 10 3 
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Fig. 1. Log- log plot of the reduced storage and loss moduli 
G'u/Cb (solid lines) and Gu/Cb (dashed lines) versus the re- 
duced frequency ujte as given respectively by the real and 
imaginary parts of Eq. 1161 (a) for T n /rE ~ 10 3 and (b) 
for Tn/i~E = 5. For the purpose of illustration we have set, 
by and large arbitrarily, vb/{CbTe) = 1.8, Di/Cs = 0.0135, 
D 2 /Di = 10 and A = -1. 



With this assumption, a reasonable phenomenological form 
for the complex modulus GT, (u>) that incorporates the ef- 
fects of both Rouse and director modes is 

GUu) = C 5 + + ^-—D 1 ) h R (-wr E ) 



TE 



2 TE 



El (El) [1 - HDl/D 2 )(T n /T E ) h R (-iLJTE) 

2 \D, 



1 + (tu/te) h R (-iu)TE) 



(16) 



G* ± (lu), on the other hand, has the same form as for an 
isotropic system, cf. Eq. J3J). The storage and loss mod- 
uli Gn(w) and G'/<(ui) implied in Eq. Ijl6|l are plotted in 
Fig. m If r n < te the frequency dependence of the moduli 
is dominated by the elastomer modes with G'^(uj) = 

and G\\(lu) = vfuj with vf = u 5 + [A 2 /(2 J T)](1 - ri/r 2 ) 2 
for lute <C 1 and G'^{tu) ~ Gj|(w) ~ y/uj for ujte ^S> 1- 
If T n 3> te, however, the structure is much richer. In ad- 
dition to the behavior for r„ < te, the storage and loss 
moduli develop respectively a noticeable plateau and dip 
at intermediate frequencies. 

Two experimental groups, the Cambridge group |14l 
ITT)irr7] and the Martinoty group [TBHTH] , have reported ex- 
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tensive measurements of GT, and G* ± on a variety of sam- 
ples. The paper of Martinoty group shows curves of G',, (w) 
and G',!(uj) that are very similar to those in Fig. ^b) for 
Tu/te = 5. Remarkably, the experimental curves are also 
compatible Fig. [TJa) where t„/t£ <C 1 if one ignores the 
low-frequency part of the theoretical curves. 

In their simplest form in which r is independent of fre- 
quency, the NED equations treat the non-hydrodynamic 
director modes but no other non-hydrodynamic modes. 
Ignoring other non-hydrodynamic modes is only justified 
if there is a clear separation of time scales with ti > Tg. 
Martinoty et al. argue convincingly that their data is in- 
compatible with this assumption. Indeed the Cambridge 
group does not disagree with this observation at least in 
some materials since in Ref. [TJ\, they indicate explic- 
itly that there is insufficient separation of time (actually 
frequency) scales for the intermediate plateau shown in 
Fig.^a) to be observed. Our view is that the experimental 
evidence is that t„ is of order te in current experimen- 
tal systems, though current data does not rule out the 
possibility of r n being as much as 5 times greater than 
te as a comparison of Fig. ^b) and the Martinoty data 
shows. Though it seems unlikely that r n ~^> te, one cannot 
even rule out this possibility entirely on the basis of the 
current experimental data, because the frequencies in the 
experiments might not be low enough to probe t" 1 if t„ 
is greater than several hundred seconds. As a comparison 
of the experimental curves and Fig.^a) shows, there is a 
danger of overlooking the intermediate plateau in G'n and 

the dip in Gjj if frequencies less than r" 1 are not probed. 

The measurements by the Cambridge group and the 
Martinoty group include 2 samples each prepared in the 
same way by both groups. For these samples, the curves 
for G' ± as a function of temperature T at different w are 
very similar. At each frequency, G' ± shows a general up- 
ward trend with decreasing T as the glass transition is 
approached. At the lowest frequency, there is a measur- 
able dip in the vicinity of the IN transition, though that 
reported by the Martinoty group is smaller than that re- 
ported by the Cambridge group. The magnitude of this 
dip decreases with increasing frequency. Both groups also 
report a decrease in G',JG'j_ at the IN transition. 

To an impartial reader, the data obtained by the two 
groups are qualitatively if not quantitatively similar. The 
two groups, however, provide different interpretations of 
their data. The Cambridge group interprets the dip in 
Gj| / G' ± as a manifestation of dynamic softening associated 
with the transition from uniaxial paranematic phase to 
the semi-soft nematic phase whereas the Martinoty group 
says that the dip is simply explained in terms of the de 
Gennes model, in which the reduction of G\i arises from 
the D\jD\ term, without recourse to the concept of semi- 
soft elasticity. 

The de Gennes theory for nematic elastomers defined 
by the total elastic energy Ti provides a correct description 
of any uniaxial elastomer regardless of how it is formed. A 
soft uniaxial elastomer formed via spontaneous symmetry 
breaking from an isotropic elastomer is merely one with 



= C 5 - Dj/(2D 1 ) = 0. Since the high-tcmperature 
phase of single-domain nematic elastomers are at least 
weakly uniaxial, it is true that all single-domain phases 
can be described by the de Gennes model without invok- 
ing the concepts of soft or semi-soft elasticity. This ap- 
proach, however, misses the fact that the vanishing of 
and the resultant soft elasticity in nematic phases formed 
spontaneously from isotropic elastomers is a consequence 
of symmetry and not of arbitrary juggling of parameters. 
A transition from a weakly uniaxial paranematic phase to 
a more strongly uniaxial nematic phase will be accompa- 
nied by a decrease in that increases as the anisotropy 
of the paranematic phase decreases. In the resultant semi- 
soft nematic phase, C§ vanishes with the anisotropy of the 
paranematic phase. Our preference is to use the concepts 
of soft and semi-soft elasticity to describe the reduction 
of that arise from the paranematic-to-nematic transi- 
tion. We agree, however, that these concepts become less 
meaningful in samples that are really above the mechan- 
ical critical point |23| at which this transition becomes 
continuous. Being above this critical point is, however, 
incompatible with the large plateau in the semi-soft non- 
linear stress strain curve [21] ■ The Martinoty data does 
seem to be consistent with being above the mechanical 
critical point. If it is, them we would like to see exper- 
iments done on materials that are softer in the nematic 
phase if they really do exist. 

The dynamics of nematic elastomers are clearly very 
complicated. They have all of the complexity of rubbers 
and gels, which have localized gel and other modes that 
we have not even discussed and eventual transitions to 
glassy phases, in addition to that of nematics. In spite of 
these complexities, the experimental results for the com- 
plex storage modulus of a number of samples of nematic 
elastomers measured by two groups are remarkably sim- 
ilar. Neither group finds a truly soft nematic phase with 
a shear modulus an order of magnitude or more smaller 
than it is in the high-temperature nearly isotropic phase, 
but their data can be described in the context of semi-soft 
elasticity. 

The above discussion summarizes what current theo- 
ries can say about the equilibrium rheological properties 
of homogeneous nematic elastomers. The only rigorous 
statement is that semi-soft nematic elastomers are uni- 
axial solids with frequency-independent shear moduli and 
viscosities in the low- frequency limit. At higher frequen- 
cies, contributions from non-hydrodynamic modes, includ- 
ing director relaxation, contribute to the complex mod- 
uli. Though experiments seem to be converging to repro- 
ducible results, it is clear that there is much more to do 
before we have a complete understanding of the dynamics 
of nematic elastomers. 

We gratefully acknowledge support by the National 
Science Foundation under grant DMR 0404670. 
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